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Abstract. We consider the approximate control of solitons in generalized Korteweg-de Vries 
equations. By introducing a suitable internal bilinear control on the equation, we prove that 
any soliton is locally null controllable, and moreover, any soliton can be accelerated to any 
particular positive velocity, after a suitable large amount of time. Precise estimates on the error 
terms and the rate of decay in the approximate null controllability result are also given. Our 
method introduces a new insight on the control of nonlinear objects, from the point of view of 
interaction and collision problems for nonlinear dispersive equations, recently developed by Y. 
Martel and F. Merle |21II22| . It can be applied, in principle, to several other models with soliton 
solutions. 



1. Introduction 

In this paper we consider the problem of controlling solitons of subcritical, generalized Korteweg- 
de Vries equations (gKdV). More precisely, we look for an internal control f = f(t,x) applied with 
the purpose of describing the dynamics of generalized solitons in the equation 

ut + (««x + «")» = /, P = 2,3,4. (1.1) 

Here u = u(t,x) is a real- valued function, and (t,x) <E R 2 . When p = 2 and / = (|1.1[) is the 
well-known integrable Korteweg-de Vries equation (KdV). 

Additionally, if only / = 0, equation is a standard gKdV equation. It has special solitary 
wave solutions called solitons, of the form 

u(t,x) = Q c {x-ct), Q c (s) := c^^Q^s), c> 0, (1.2) 

with 

0(s) ■ P + 1 

U " [2cosh 2 (i(p- l)s) 

The parameter c > is usually denoted as the scaling, or in a equivalent way, as the velocity of 
the soliton. Inserting the previous profile in (|1.1[) (recall that / = 0), one has that Q c > satisfies 
the nonlinear ODE 

Q"-cQc + Q p c=0, Q c eH l (R). (1.3) 

Moreover, standard conservation laws for (jl.ip at the ii" 1 -level are the mass 

M[u](t) ■=- [ u 2 {t,x)dx = M[u](0), (1.4) 

and energy 



E[u](t):=- [ ul(t,x)dx — [ u p+1 (t,x)dx = E[u](0). (1.5) 

2 Jr P + 1 Jr 

A satisfactory Cauchy theory is also present at the H 1 level of regularity, see e.g. Kenig-Ponce- 

Vega [12]. The condition on p is necessary to get global existence for general iJ 1 -data, cf. the 

paper by Martel and Merle [19] for the critical case p = 5. 
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The control problem for the non inviscid KdV equation in a finite length interval has been 
extensively studied in the last twenty years, starting from the works of Zhang [39] . Russell and 
Zhang in [351 134) , and [36j for a system with periodic boundary conditions and with an inter- 
nal control. In the case of a boundary control, see [36] and [37]. Concerning the non periodic 
framework, Rosier studied [31] the controllability of the KdV equation posed on a finite interval 
of (0, L), under homogeneous Dirichlet boundary conditions and a control acting on the Neumann 
data at the right end-point of the interval. In particular, Rosier shows that if the length L does not 
belong to a set of critical values, both the associated linear and the nonlinear systems are exactly 
controllable. When L is critical, the linear system is not controllable because of the existence of 
a finite-dimensional subspace of unreachable states. In this case, the exact controllability of the 
KdV equation, in the case of critical domains, has been proven by Coron-Crcpeau [B], Cerpa [7], 
and Ccrpa-Crepeau [H]. Concerning the exact boundary control problem in the half-line, see e.g. 
the work of Rosier [35] . 

In this paper, unlike the previous results, we are interested in the study of a control problem 
associated to a given gKdV soliton posed in the real line. Let us explain in more detail the 
problem. Given an initial datum u(to, x) = uq(x) = Q c (x — cto — xq) of soliton type, our objective 
is to introduce a control / in the gKdV equation during an interval of time [0, T], with the 

purpose of accelerating the soliton to a new soliton state, with a different (positive) velocity. With 
no loss of generality, we can assume to = xq = and the initial velocity satisfies c = 1. In other 
words, our goal is to determine sufficient conditions on / to ensure that, given any final scaling 
Cf > 0, the system with initial datum uq evolves to a soliton of the form Q Cf , up to some 
small error terms, in a suitable time of interaction T > 0. Moreover, we also want to estimate the 
position of the soliton, compared with the theoretically expected position ~ CfT '. 

We will assume that the interior control / is given by the bilinear control (or feedback law) 

f(t,x) = a(t,x)u(t,x), 

with a an internal potential satisfying the a priori assumptions 



In other words, our control will use some explicit information of the soliton solution at each time, 
such as the scaling and position parameters. These problems have been also considered in a more 
physical context by Kaup- Newell [11], Grimshaw [10], Ko-Kuehl [13], and Lochak [15] . 

Therefore, in what follows, we consider the initial value problem 



where a is an unknown control. Our first result states that any gKdV soliton is locally null- 
controllability for sufficiently large time. 

Theorem 1.1. Any gKdV soliton is approximate null- controllable in large time. More precisely, 
given any <5, <5o > small, there exists Eo(8,8o) > such that, for all < e < eo, the following 
holds. There exist a time T = T £ > 0, and a smooth control a = a e (t, ■) G i°°(M) n L 2 (R), defined 
in [0,T] , such that the unique solution u = u e (t) of in C([0, T], iJ 1 (R)), with initial condition 
Uq(x) = Q(x), satisfies 



Finally, one has T ~ e 1 . 

As far as we know, this is the first (partial) result on controllability of solitons in unbounded 
domains, where dispersion plays a key role in the dynamics. Previous results are related to the 
study of the ground state of the linear and the nonlinear problem in a finite interval, see e.g. the 
works of Lange and Teismann [14] , Beauchard and Mirrahimi [2] [25] , or the control of a quantum 
particle under the action of a well shaped potential, obeying the linear Schrodinger equation in 
a bounded interval (Beauchard-Coron [3]), and Crepeau [9] in the KdV case. In this paper we 



a{t, •) G C* 3 (M) n L 2 (R) n i°°(R). 



(1.6) 




(1.7) 



\\u(T)\\ mm <5. 
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study a nonlinear object instead; localized solitons are present due the non compact character of 
the domain, and the strength of the nonlinearity. 

It turns out that Theorem 11.11 is consequence of the following deeper result, a large time ap- 
proximate controllability of the initial soliton Q of scaling one, to any final scaling c/ > 0, c/ 7^ 1 
(the case c/ ~ 1 is trivial). As previously stated, we pick any do > small. 

Theorem 1.2. Let c/ > 0, c/ 7^ 1. There exists eo(^o,c/) > such that, for all < e < £0 the 
following holds. There exist a time T — T e > 0, a smooth control a — a £ (t, ■) £ i°°(R) fl L 2 (W) 
and a smooth translation parameter pit), both defined in [0, T], and such that the unique solution 
u = u £ (t) of tl.l\) in C([0, T], H 1 (M.)), with initial condition Uo(x) = Q(x), satisfies 

\\u{T) - Q Cf (- - p(T))\\ HHK) + \p'{T) - Cf \ < V~e. (1.8) 

Proof of Theorem 11.71 assuming Theorem \1.2\ The proof of Theorem 11.11 follows from Theorem 

4(p-l) 

11.21 and the subcritical character of solitons for p < 5. Indeed, just take any < c/ < S 5 -p 
and So < <5 in Theorem 11.21 Note that Theorem 11.11 holds even without destroying the soliton 
structure. □ 

Some comments about Theorem II. 21 

Remark 1.1. First of all, we point out that the control is not compactly supported, but it satisfies 
the following properties (see Section [2] and (|2.6p for more details): 

(1) It is exponentially decreasing in any moving region far away form the soliton (in other 
words, it moves with the soliton); 

(2) It has slow variation in space, which actually explains the large time needed to drive the 
dynamics. 

Therefore from (|1.8[) . the solution u(t) satisfies, after some direct verifications, the generalized 
observability inequality 

WQWlhr) < I I \u{t,x)\ 2 dxdt. 
Jo Js. 

Remark 1.2. Second, as for the final position and time of control T are concerned, we obtain 
estimates of the following orders: for any 80 > small but fixed, 

T ~ e- 1 - 50 , \p{T) - c f T\ < £ - 1/2 - So , 

although the relative error in the last estimate is (^(e 1 / 2 ). The relative weakness in e of the 
last estimates and the bound (|1.8[) is mainly due to the emergence of dispersive tails behind the 
soliton solution as the control acts; this phenomenon has been observed in several interaction 
problems involving gKdV equations, starting from the formal arguments in (TTJ [TB] HD] , and the 
more rigorous treatment given in |211 l2"2l l2l)l |2"7] . This phenomenon does not appear in the case of 
nonlinear Schrodinger equations, where one expects better estimates (see e.g. [27]). Heuristically 
speaking, the lack of control on the position, compared to the extremely accurate control on the 
velocity, could be associated to a form of uncertainty principle for solitons, regarded this time as 
approximate point particles. 

Remark 1.3. A necessary condition to obtain an estimate as in (|1.8|) is the lack of conserved 
quantities (see Proposition ^. 1[) . Indeed, it is not difficult to see that the soliton at time T has lost 
or gained, depending on the sign of a and c/, a nontrivial O(l) amount of mass (|1.4j) . A similar 
study can be applied to the case of the energy (|1.5[) . with similar conclusions. It is important to 
stress that, since solitons are stable under small H 1 perturbations [U 21 135] and the equation is 
not integrablc unless p = 2 or p = 3, it is expected that the above result only holds if we introduce 
a sufficiently slow varying potential. 

On the other hand, problem can be also regarded as a stabilization problem. In that sense, 
the recent literature concerns with the decay of solutions posed in a bounded interval [33 , 23 , 5 , 15 , 
or the half line by Linares and Pazoto [TTJ [TSJ [5S] , and numerical schemes for the critical case p — 5 
(Pazoto ct. al. [H]). As for the decreasing mass case, and the approximate null controllability 
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result stated in Theorem 11.11 we have the following additional local stabilization result, without 
destroying the soliton: 

Corollary 1.3. Under the assumptions of Theorems 1 1. 1\ and \1.2l there exist C, po > ; indepen- 
dent of e € (0,£o) such that, for all t G [0,T] 

IM*)||j?i(H) <C(^ + e-^ et )\\Q\\m( R) - (1.9) 

1.1. About the proofs. Our proofs do not involve the usual methods employed in control theory, 
requiring the study of the linear problem, unique continuation properties and/or Carlcman esti- 
mates. In order to study genuine nonlinear objects such as solitons, nonlinear dispersive methods 
are somehow needed. In particular, a suitable global wcll-poscdncss theory in the energy space 
for solutions of (|l.lj) in the real line requires modifications on the arguments of the fundamental 
work by Kenig, Ponce and Vega (T2], in order to deal with the unbounded domain case. Second, 
our control is explicitly constructed, with the following properties: (z) it has a slowly varying 
character, determined by the parameter e; (ii) it is localized in a moving region of size O(l), and 
is of strength 0(e) (but it is not compactly supported), and (Hi) the corresponding slowly varying 
part induces on the soliton parameters a finite dimensional dynamical system which governs the 
whole dynamics. Concerning the time of control, since the dynamics is slowly varying, the time 
of interaction is 0(e~ l ) at least; a control introduced in a smaller window of time could destroy 
the soliton. 

The second step of the proof is the following: since the introduction of the control induces on 
the soliton the action of an external potential, we can think such an interaction as a slowly varying 
collision between both objects. In Section |3l we construct an explicit approximate solution which 
describes the interaction, up to certain order of accuracy in e. This solution u(t) has the form 

u(t,x) ~ Q c(t) {x - p{t)) + eA(t,x- p(t)), 

where e > is a small, artificially introduced parameter, and (c(t), p(t)) are suitable scaling and 
translation parameters, depending on time. The parameters follow a suitable approximate finite- 
dimensional, slowly varying in time dynamics, determined by the action of the control, described 
as follows: 

c'{t)~a(t,ep(t)), p'(t)~c(t)+0(e), 
where a is the control introduced in (|2.6I) . We choose carefully a such that the evolution of this 
system leads to the desired final velocity, at time T ~ e -1- * 50 , 

c(r)~ C/ , p(T)~c f T + o(T), (1.10) 

however, a better control on the position has escaped to us. 

Concerning the function A, it corresponds to a first order correction term with support of size 
0(e~ 1 ) in the variable x — p(t), and _L°°-norm of order O(e~ loe ^ p ^), for some constant 70 > 0. 
Therefore, A is a phantom term that disappears after the interaction, but which allows to improve 
the accuracy of the approximate solution. Additionally, A is in principle only bounded, but not 
localized^] therefore we introduce a suitable cut-off function to recover a finite mass solution. The 
error associated to this approximation is measured in terms of the L^°H^. norm, and it has to be 
small enough in order to take into account the large time of interaction. In our case, we are able 
to prove that during the whole interaction, one has 

remainder - £ 3 / 2 e - 7o,r|p(t)l , f g 

(see (|3.45[) ). therefore the propagation of this error during a time interval T formally leads to 
the bound 0(y/e) in Theorem 11.21 We remark that this method has been recently applied, in a 
different context, to several interaction problem, notably the two-soliton collision by Y. Martel 
and F. Merle (2UE2], an d the interaction of solitons with a potential [26], [27]. See also (30] for a 
related soliton-potential problem, in the easier case of the cubic nonlinearity, and for which the 
term A is not needed. 



In principle, A models a dispersive tail behind the soliton solution. 
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The third step of the proof is the following. In order to control the dynamics of the error 
terms, we introduce a suitable Lyapunov functional (Section 2]), adapted this time to the genuine 
nonlinear dynamics of the problem (see e.g. [21j). This functional has very small variation in 
time, provided we control the size of some time dependent parameters of the soliton solution. We 
avoid that problem by using sharp virial estimates, in the spirit of |20j . After this point, we can 
close the main argument by proving rigorously that the error terms can be assured to be smaller 
than 0(y/e), during the whole interaction region. 

The final step of the proof is a rigorous analysis of the time parameters (c(i), p(t)) of the soliton 
solution, in order to recover (jl.lOp . We prove that at time t = T, the solution has the desired 
behavior, up to an error of 0(\/e), concluding the proof of Theorem 11.21 Finally, the proof of 
Corollary 11.31 follows after a detailed study of the scaling parameter c(t). 

The weakness of our approach is precisely in the approximate character of the controllability 
property, and the large time needed to reach an approximate final state. We believe that our 
results can be improved by adapting, to this case, the standard and complex machinery of control 
theory. Moreover, we believe that the moving profile of the support can be chosen to be compactly 
supported. However, we do not expect an exact controllability result, in view of the following 
remark. 

Remark 1.4. Note that if (|1.7[) is exactly controllable to zero in finite time, say u(T) = for some 
bounded, smooth control a(t, x), then using the reversibility in time of the equation, we have that 
v(t, x) := u(T — t, —x) satisfies a slightly different equation, 

«t + (vxx + v p ) x = b{t, x)v, «(0) = 0, 

for the potential b(t, x) := a(T — t,—x). It turns out that, under standard assumptions on the 
solvability of the Cauchy problem associated to v, the unique solution to the above problem is the 
identically zero solution, a contradiction. We thank Sylvain Ervedoza for this remark. 

We point out that in order to describe the dynamics in a time of order O(l), one needs formally 
a large control; in particular, it should be unbounded in space (more precisely, linearly growing in 
space). However, even the local in time Cauchy problem for such perturbations becomes a very 
difficult problem. 

Finally, some words about the organization of this paper. In Section^ we introduce the explicit 
control system, the finite dimensional dynamical system and the corresponding local and global 
well-posedness theory. In Section [3] we construct an approximate solution to a given order of 
accuracy. We continue this process up to the moment when we find an infinite mass correction 
term, which is up to date the best mathematical description of the dispersive tail behind the 
soliton, originated by the application of the control. Section |4] is devoted to the introduction of a 
Lyapunov function, modulation theory and a key virial identity in order to control the dynamics 
of the oscillatory terms. Finally, in Section [6] we prove the main theorem. 

Notation. Along this paper we use the convention A < B if and only if there exists K > 0, 
independent of e, such that A < KB. Additionally, 7 and K* will denote special positive constants, 
still independent of e, to be worried about. Finally, <S(R) denotes the Schwartz's class on R. 

Acknowledgments. The author would like to thank Gunther Uhlmann and Axel Osses for 
their kind invitation to the PASI-CIPPDE 2012, Inverse Problems and PDE Control, held in 
Santiago-Chile, and where this project was originally conceived. The author is very grateful of 
Eduardo Cerpa, for many useful comments and suggestions to a first draft of this paper. 
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aoc :=-^log C/l p = 2, aoo := } P ^ (l-cf 1 ), p = 3,4. (2.2) 
A2 A p (p — I) 

We introduce now the control a(t,x). Given any e > small, we consider a smooth function ao 
satisfying the following properties 

'a E ^(IjflL^fR), 

|ao(a?)| < e 70 ", for x < -1, | aoo - a (x)\ < e~^ x for x > 1, 
14^(^)1 ^ e-T°N, i£R, fc = l,2,3, 
,a' (i) > if > 0, Oq(x) < if < 0, 

for a fixed, positive constant 70. Note that with our choice, 

IklU = ko|. (2.4) 

Let < c m :— imin{cj,l}, c M := 2max{l,c/}, and (c (t), p (t)) E R + x R be a set of C 1 
parameters defined in some interval of {t > 0}, with the following uniform, a-priori constraints 

< c m < co(t), p' (t) < c M , 

I^WI + ^W-coWl^e 1 / 2 " 50 , (2.5) 
JcoW-ll + lpoWl^e 1 / 2 - 5 ", 

for some <5o > small. Finally, define 

a(t,x) := -ea' Q (ex)Q C0 ( t) (x - po{t)), (2.6) 

where Q c is the solution of (|1.3[) . Note that a, as a function of the pair (co,po), is a contraction 
in any small vicinity of that point. Indeed, we have 

|o(co,Po) - o(ci,pi)| < e(|c - ci\ + \p Q - px\). 

Let us remark that this control will take into account some important information of the soliton 
itself, namely the scaling c(i) and the position p(t), and it is in some sense of nonlinear character. 
In terms of applications, these two parameters can be easily described by a suitable approximate 
dynamics, described by (|4.2[) - (|4.3p . 



It is not difficult to check that this control satisfies the following space-time bounds 

ll a IU°°([0,oo)xR) + ||a:r:r||L=°([0,oo)xR) ^ £■ 

Under these estimates, we claim that the Cauchy problem associated to (jl.ip is locally well-posed 
in a subspace of iJ^R). 

Proposition 2.1. Under the assumptions \2. 3\) , \2. 5\) and &2.6\) , the initial value problem 

is locally well-posed in H 1 (R). Moreover, the mass M[u](t) and energy E[u](t) defined in \H ) 

and il.5\) satisfy the relations 

d t M[u](t) = - I a(t,x)u 2 , d t E[u]{t) = \ [ a xx u 2 + / au p+1 - [ au\. (2.7) 
Js. 2 J K J M J R 

Finally, the map (c, p) — > u(c, p) is a continuous contraction in a small neighborhood of (co, po). 

Remark 2.1. Later we will prove that our solution is well-defined, for all t < T ~ e^ 1 ^ 6 ", as a 
consequence of the stability property (fTT 



Proof. This result is classical, see e.g. Bona-Zhang [S] for the KdV case, or Merle- Vega [23] in 
the mKdV case. For the sake of completeness, we sketch the main details. We use the machinery 
developed by Kenig, Ponce and Vega [12] in order to prove local well-posedness for gKdV in low 
regularity Sobolev spaces. Since we only need an H 1 local theory, our proofs will be simpler than 
the original ones. 

Recall that we want to solve 

ut + u xxx = -[pu p ~ l u x + a(t,x)u], u(0) = uq e ff 1 (R). 
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If we denote by e * the free Airy propagator, we have to solve fixed point problem 

ft 

u(t) = T[u](t) := e- ta *u Q - / e^*" 5 ^ [p^" 1 ^ + au](s)ds := To[u ](t) + Ti[u}(t). 



n 



Note that, since we have chosen (c, p) following (|2.5|) . 

/ sup a 2 < e 2 / a,Q(ex)dx < e, sup ||a||z,«.( R ) < s. 

JlO<t<l JR 0<t<l 

Using [12j Theorem 3.5] and the maximal function estimate |T2l (3.9)], we have, for any S £ (0, 1), 

sup ||da;7iM||L2( E ) < / Wpu^Ux + au|| L 2 (0 < t < 5) cfa; 

0<t<S 

<( f sup \u\ 2p ~ 2 dx) 1 |ju a; |jL2 ([0 .5 ]xK) + ( ( sup a 2 dx] 1 ||u|| L 2 ([0iS]xR 

V JlO<KS ' ' V JlO<KS J 

< S 1/2 \\u\\ p L ^ ao syHim +v / e-S 1/2 ||u|| LO o(( 0iS );L2( ffi)) . 



On the other hand, 



|7iMIU 2 (iy < / [||" p 1 Wx||l2 (Kx) + ||au|| L 2 (Rx) ]di 

" S - 1/2 



<5 sup (/ \\u\\ p L J (R )) |K|| £ 2 (Ric ) +S sup ||a|| L ^ (Rx) ||w|| L 2 ( 



0<t<S y Jo v ' 0<t<S 

It is not difficult to check that these estimates give that, for S small, T a ball of H 1 into itself. 
The contraction follows in a similar way. 

Let To > be the maximal time of existence of a solution u(t). It is not difficult to check that 
the mass and energy (|1.4|) - ()1.5j) satisfy, for all t <G [0,T ), 



d t M[u](t) = - / a(t,x)u 2 < sM[u](t), 

therefore M[u](t) < e Cet . On the other hand, the energy (|1.5[) satisfies the relation 

d t E[u](t) = - I a xx (t,x)u 2 + / a(t,x)u p+1 ~ / a(t,x)u 2 x . 
2 Jr Jr Jr 

It is important to stress that in the case a > (c/ < 1) and p = 3, this estimate gives global 
existence for u. However, in a general case, a better estimate is somehow needed. □ 



3. Approximate solution 

Given any e > and Sq > small, we introduce the time of interaction 

T^minlTo,^ 1 - 50 }, (3.1) 

where To > is the maximal time of existence of the solution u(t) with initial condition Q(x). 

Let t G [0,T]. In what follows, we fix a couple of dynamical parameters (c(t),p(t)) satisfying 
(|2.5[) in some subinterval of {t > 0}. Additionally, we will need the extra assumptions 

\c{t)~c {t)\ + \p{t)-p {t)\<e l ' 2 - s \ (3.2) 

valid for all t e [0,T]. 

Now we define the modulated soliton solution as follows 

y:=x-p(t), R(t,x) := Q c {t){y)- 
Finally, we introduce the approximate solution 

u(t, x) := R{t, x) + w(t, x), w(t, x) := ed(t)A c(t) (t, y), d(t) := a' (ep(t)); (3.3) 
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for some L°°(R) function A c (t, ■), to be found later. In order to simplify some computations, we 
will assume that for c m < c < cm and t fixed, A c (t, ■) satisfies the estimates 

A' c {t, •) G L 2 (R), d c A c (t, ■) G L°°(R), (3.4) 

that will be verified later. Finally, we define the scaling operator 

AQc(y) := 9 C ,Q C , I (i/) = - \-L.Q c (y) + \yQ' c { y )\ G 5(M). (3.5) 
lc - c c Lp — 1 2 J 

We want to measure the size of the error induced by inserting u as defined in (|3.3[) in the 
equation (fl~7)) - (j2~6l) . Let 

S[u](t,x) := ii t + (u xx + u p ) x + sa' {ex)Q Co ( t )(x - p (t))u. (3.6) 

From (|3.2|) we have 

Qc (t)(a: - po(t)) = Q c (t)(x - p(t)) + &Q c (t)(x - p(t))(c(t) - co(t)) 

+Q' c(t) {x - P (t))(p(t) - p„(t)) + Oh^r)^ 1 " 250 ). (3.7) 

Our first result is the following 

Proposition 3.1. Let (c,p) be satisfying A2.5\) . There exists a function A c G L°°(R) smc/i £/iai u, 
defined in \3. 3]) . satisfied 

S[u](t, x) = (c'(t) - eh{t) - 0(e 3 / 2 - s "))d c u-(p'(t) - c(t) - ef 2 (t) - 0{e^- s °))d v u + S[u], (3.8) 
where fi(t) = /i(c(t), p{t)) and /2(f) = /2(c(i), p{t)) are given as follows 

p 2(5-2p) 

fi{c,p) := -\ p a {£p)cp-i , f 2 (c,p) := p p a {ep)c , (3.9) 
/or some G R, m'tt ^3 = 0. Moreover, 

\\~Sm)\\m { y>-i)<^e-^W)\ +£ \ (3.10) 

and 



Q c S[u] 



yQ c S[u) 



< e 2 e - e7 |p«l + £ 3 j (3 n) 



Proof. We follow the strategy described in [27]. Suppose that the parameters (c(t),p(t)) satisfy 
(j2~5l) and ([321). From we have 



5[u] = I + II + HI, (3.12) 
where (wc omit the dependence on t and x) 

l:=S[R], 11 = 1%) :=w t + (w xx +pR p ~ 1 w) x +ea' (sx)Q Ca (x- p a )w, (3.13) 
and for p = 2, 3 or 4, 

III:= {(R + w) p - R p -pR p ~ 1 w} x . (3.14) 

Recall that w is given by (|3.3p . In the next results, we expand the terms in p,12|) . Note that 
R(t,x) = Q c (t) (y) and y = x - p(i). 

Lemma 3.2. 

I = F*(t,y)+eFl(t,y)+e 2 Fl(t,y), 

where 

Fl{t,y) := (c'(t) - eh{t))d c R{t) - (p'(t) - c(t) - ef 2 (t))d x R(t), 
fi(t) and /2(f) are given by LS.9\) , and 

F*(t;y) ~h(t)AQ c (y) + a' (ep)Q 2 c (y)-f 2 (t)Q' c (y)- (3-15) 
Finally, for all t G [0,T], one has \\F*(t, -)\\h^(s) £ e'^P^ + e. 



2 The first two terms in 113.811 are often referred as the finite-dimensional dynamical system associated to the 
soliton dynamics. 
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Proof of Lemma \3.2\ We have 

I = Rt + (R xx + R p ) x + ed^{ex)R 2 = c'AQ c - p'Q' c + + (QP)' + ea' (ex)Q 2 c . 
On the other hand, note that via a Taylor expansion, 

a' (ex)Q 2 c = a' (sp)Q 2 c + sa'^ep)yQ 2 c + O m{M) (e 2 ). 
Therefore, using the equation satisfied by Q c , namely, Q" — cQ c + Q p = 0, one has 

I = c'AQ c - (pf - c)Q' c + - cQ' c + (QP)' + ea' Q 2 c + e 2 a^yQ 2 c + HHm (e 3 ) 
= (c' - e/i)AQ c ~(p'-c- ef 2 )Q' c + e[/iAQ c + a' Q 2 c - f 2 Q' c ] + e 2 F 1 c {t 1 y), 
with F} := a'^yQl + O h i {r) {s). Moreover F*(t, •) £ 5(R) and \\F*(t, -)\\m(R) < e~ s ^P^ + e. □ 
Now wc recall the linearized elliptic gKdV operator. Fix c > 0, p — 2, 3 or 4, and let 

Cw := -w yy + cw - pQ p c r 1 {y)w, where Q c (y) '■= c~Q(^/cy). (3.16) 
Here w = w{y). 

Lemma 3.3. Suppose that A c satisfies \3.J$ . Let w given by II 3.3)) . Then the following expansion 
holds: 

II = (c' - eh)d c w -{p 1 -c- ef 2 )w y - {Cw) y + e 2 [a%cA c + ha' d c A c ] + e 2 F*\t; y), 

with 

F}\t; ■) = ea'J(ep(t))f 2 (t)A c + a'^ep(t))(p> (t) - c(t) - ef 2 {t))A c + O^^-^"^), 
for some fixed 7 > 0. 

Proof. Let D := D c (t, y), y = x — p(t), be a general, smooth function. We compute 
11(D) := D t + (D xx + pRP^D)* + ea' {ex)RD. 

We have 

11(D) = c'(t)d c D + D t - (p'{t) - c(t))D y + [D yy - c{t) D + pQ p c - 1 D] x + ea , (ex)RD 

= D t - (CD)y + (c'(t) - eh{t))d c D - (p'(t) - c(t) - ef 2 (t))D y + O(ea' (ex)Q c D) 
+ef 1 (t)d c D-ef 2 (t)D v . 
We apply this last identity to the function w = ed(t)A c (y). We have 

II(w) = ed'{t)A c - ed(t)(CA c y + (c'(t) - eh(t))ed{t)d c A c - (p'(t) - c(t) - ef 2 (t))ed(t)A' c 
+e 2 d(t)f 1 (t)d c A c + HHR) (e 2 e- e ^\) 
= (c'(t) - eh{t))ed{t)d c A c - (p'(t) - c(t) - sf 2 (t))ed(t)A' c - ed(t)(CA c )' 

+e 2 [e- l d'{t)A c + d(t)h{t)d c A c ] + H i m (e 2 e-^"^). (3.17) 
(Recall that A' c £ S.) Now we use the fact that d(t) = a' (ep(i)) to compute d'(t). We have 

d'(t) = e4(e P (t))p'{t) = ea'^ep{t))c{t) + e 2 a^ep(t))f 2 (t) + ea^(ep(t))(p' (t) - c(t) - ef 2 (t)). 
Replacing in (|3 . 1 T[) we conclude. □ 
Lemma 3.4. Suppose that A c satisfy \3.$ . Then 

\U = O mm {s 2 e-^\). (3.18) 

Proof. First of all, define III := (R + w) p - R p - pR p ~ 1 w. Then, 

'e 2 d 2 (t)A 2 c if jj = 2; 

III = < e 2 d 2 (t)A 2 c [3Q C + ed(t)A c ] if p = 3; 

e 2 d 2 (t)A 2 c [6Q 2 c + 4ed(t)Q c A c + e 2 d 2 (t)A 2 } in the case p = 4. 

Thus taking space derivative we obtain (I3.18|) (note that (A p )' £ S because A c satisfies (|3.4[1 ). □ 



in 
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Now we collect the estimates from Lemmas l3.2[ 13.31 and l3. 41 We obtain that, for all t in a given 
interval, 

S[u] = (c'(t) - eh{t))d c u - (p'(t) - c(t) - ef 2 (t))u y + S[u], (3.19) 

with 

S[u] = e[Fi(t, y) - d{t){CA c ) y ] (3.20) 
+ e 2 [a^(ep(t))cA c + f ia ' (ep(t))d c A c ] + s 2 a%{ep(t)){p'(t) - c(t) - ef 2 {t))A c (3.21) 
+ e 3 ai>(ep(t))f 2 (t)A c +e 2 Hlm (e~^r^+e). (3.22) 
In addition, fi(t),f 2 (t) are given by (|3.9|) . and 

Fi := Fl = h(t)AQ c + a' {ep{t))Q 2 c - f 2 (t)Q' c , (3.23) 
(cf. (|3.I5[1 ). Moreover, it is not difficult to check that, for any t, one has 

F 1 (t,y)Q c (y)dy = 0. (3.24) 

The next step is the resolution of the linear differential equation involving the first order terms 
in e. Indeed, from (|3.20p . we want to solve 

d{t){CA c ) y {y) = Fi(t,y), for all y € K, and t fixed; 

with d(t) given by (|33]l . Note that from (|33|) and (|3T2^|) one has 

Fx(t; y) := a' Q [-X p c^AQ c + Q 2 C ] - f 2 (t)Q' c , 

P n 2(5-2p) 

= a [-A pC — KQ C + Q 2 C -p p c^i-Q' c ] =: d(t)Ji(*,i/). (3.25) 
Therefore, we are reduced to solve the following simple problem, 

(CA c ) y ( V )=F 1 {t 1 y), (3.26) 
with Fi defined in (|3T25|) . and from ([3~24]) . 

F 1 (t,y)Q c (y)=0. 

Let us recall the following results (see e.g. |21j): 

Lemma 3.5. The operator C defined (on L 2 (M)) by p.I6p /ias domain H 2 (M), it is self-adjoint 
and satisfies the following properties: 

(1) The kernel of L is spanned by Q' c . Moreover, AQ C defined in \3.5\) satisfies C(AQ C ) = —Q c . 
Finally, the continuous spectrum of C is given by cr con t(C) = [c, +co). 

(2) For all h = h(x) polynomially growing function such that J R hQ' c = 0, there exists a 
unigue polynomially growing function h such that J R hQ' c = and Ch = h. Moreover, if 
h is even (resp. odd), then h is even (resp. odd). 

(3) For h e H 2 (M.), Ch G «S(R) implies h G «S(R). 

We introduce the following function, with the objective of describing the effect of potential on 
the solution. Let c > and 

¥>(*):= -tS^, <Pc(x):=-y?=V~c<p(V~cx)- (3.27) 
Q(x) Q c 



Note that (p is an odd function, with 



lim ip(x) = ±f ; ip w e 5(R), k > 1 . (3.28) 

x— > ±oo 

We recall the form of the solution A c that we are looking for. We seek for a bounded solution 
satisfying 

A c(t )(y) = Pc{t)(y c {y) - v^W) + Mv) + »c(t)Q' c (y) + 6 c (t)AQ c (y), (3-29) 
for some f3 c (t), p c (t), 5 c (t) G M, ip c defined in (|3.27[) . and A c G S. The parameters p c , S c will be 
chosen in order to find the unique solution A c satisfying some orthogonality conditions. 
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Lemma 3.6. Suppose (c(t),p(t)) satisfying 
unique solution A c n\ = A c ^{y) of 



id fi{t), f2(t) given by i3.9\) . There exists a 



such that, for every t, 



(£A c(t) ) v (y) = h{t,y), 

A c ( y ) ■= MMy) ~Vc) + Mv) + McQ' c + S c AQ c {y), 
limA c = -2 Vc/3 C ; \A c (y)\ < Ke~ lv , as y -> +00 , 



with A c £ S. In addition, we hav^ 

1 



Pc{t) := 
Moreover, A c satisfies 



2c 3 / 2 (i) 



F x {t,y)dy £ 0, |/3 c (t)| + |^ c (t)| + |* c (f)| < 1. 



/ A c Q c (y)dy = / A c yQ c (y)dy = 0. 



(3.30) 

(3.31) 
(3.32) 

(3.33) 
(3.34) 



Proof. First of all, the existence of a solution A c G L°° (R) of the form (|3.31l) for this equation was 
established in [26], provided 

r F 1 (t,y)Q c = 0, 



which is indeed the case (cf. p.24[) V Note that the inclusion of the term proportional to f2(t)Q' c 
in (|3.23[) induces the new term 5 C AQ C in (|3.3ip (recall that from Lemma [3.51 (CAQ C )' = —Q' c .) 
Furthermore, the limits in fl3.32[) are straightforward from (|3.28[) . 

On the other hand, we choose the terms fi c and S c in order to satisfy (|3.34[) . Since we do not 
explicitly know A c , we need another method to compute explicitly f2{t) (and therefore, 6 c (t)). 
Indeed, multiplying p.30|) by f_ AQ C € L°°(R) and integrating, one has 



(CA c ) y / " AQ C = f Fx f V AQ C . 

> JS. J -00 



(3.35) 



Integrating by parts, we get 
(CA C ) / AQ C 



(£A C )y 

+00 



AQ C 



AQ C CA C 



JL = 0. 



0. Therefore, from (|3~25)) . 

y 



Using dHU), we have (CA C ) AQ 

-/ 2 J QcAQ c = a' Q J[- X p c^AQ c + Q 2 C 
A simple computation shows that 



AQ C . 



-0f 2 c 



261-1 



1, _E_ 



AQ C 



Ql / AQ C 



(3-p) 1. (3-p) 



2(p-l 



■a n c p- 1 



~4 Ap ~^T 



QY 



Using that 8 = — |, we finally obtain 



(5-p) 
u p a' (e/a 

as desired (cf. (|3.9p ) (note that ^3 = 0) 



4 p p- 1 



QY 



Q 



( Q 2 ) 



aia-2|>) 



3 Note that /3 C = implies A c G L 2 (ffi). 
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Now, let us prove (|3.33[) . Indeed, from (|3.30p . integrating over R and using (|3.32[) . we get 

2p c cy/c=aA c (-oo)=CA c (+oo)-CA c (-oo)= / F u (3.36) 

Jr 

which gives the value of /3 C . Moreover, 



F^c 26 / Q^\ pC 2 %-^-t) / Q 



c 29 



1, 

K p-1 2' 



value: note that from the identities Q" = Q — Q 2 and Q' 2 = Q 2 — iQ 3 , one has 



for p = 3, 4. The case p — 2 requires more care, but a simple computation gives a nonzero final 

!< 

Q\ I Q 3 = | ! Q 2 = f f Q. 

JR JR JR 

Therefore ( J M Q 2 ) 2 - § / R Q / R Q 3 = (1 - | ) ( f M Q) 2 > 0. 

Let us now describe the dependence on c of the solution A c . From (|3.25[) (see also Lemma 4.5 
in [IS]), one has 

2 32p-13-llp 2 ~„ 

Fi(t,y) = cv=iFl(y/cy) + c'O-w-sp) F{(y/cy), 

where 

Fl{y) := -X P AQ(y) + Q 2 (y), F 2 (y) := Q' . 
Moreover, Claim 3 in [26] allows to conclude that A c satisfies the following decomposition: 

A c (y) = c$^Al(y/cy) + d'^-sp) Al(y/cy), (3.37) 

with A\ bounded solution of {LA\)' = F/, and A\{y/cy) G 5(M). Moreover, one has (i})' 6 5. 
Using this decomposition we have that 9 C A C has the same behavior as A c : it is bounded, it is not 
L 2 -intcgrablc, and satisfies lim +oc d c A c — 0, hm_oo d c A c ^ 0. The same result holds for d 2 A c . □ 

Having solved the linear problem, from (|3.19|) and (|3.20j) we have 

S[u](t, x) = (c'(t) - efi(t))d c u - (p'(t) - c(t) - ef 2 it))u y + S[u}(t,x), 

where S[u] given in (|3.2ip - (|3.22|> is now of second order in e. Note that from (|2.5p we have 

\p'(t)-c(t)-ef 2 (t)\<l, 

uniformly in time. Moreover, using the decay on the right of A c (see (|3.32I0 . estimate (|3 . 10[) is 
direct. In addition, from Lemma [3.61 we have Q3.4p . and fi(t) and fz(t) are well determined by 
(EH) . Finally, from ([3T2"l"]) - (pT2"2"]l one has ([3~TT|) . These facts prove Proposition 1331 □ 

The next results are similar to those proved in [55] HZ], but for the sake of completeness, we 
include them. Recall that u does not belong to L 2 (R). In order to solve this problem, consider a 
cut-off function r\ £ C°°(M) satisfying the following properties: 



(3.38) 



I < T](s) < I, < t]'(s) < 1, for any s G M; 
1 r/(s) = for s < —1, ?j(s) = 1 for s > 1, 

and define 

:=?7(ey + 2), (3.39) 
and for w = w(t,y) the first order correction constructed in Lemma 13.61 redefine 

u(t, x) := rj £ {y)u(t, x) = rj e (y)(R(t, x) + w{t 7 x)), (3.40) 

and similarly for R(t) and w(t). Note that, by definition, 

u(t, x) = for all y < ~. (3.41) 

The following Proposition deals with the error associated to this cut-off function, and the new 
approximate solution u. 
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Proposition 3.7. There exist constants e ,7 > such that for all < e < Eq the following holds. 

(1) Consider the localized function u(t) = R(t) + w(t) defined in i3. 39\) - [3.40j l, for t in a given 
interval. Then we have 

(a) i 2 -solution. w(t,-) g H 1 ^), with 

\Htr)\\ H i(m)<^e- 1£im - 

(b) Almost orthogonality. 



w{t,x)Q c (y)dx 



yw{t, x)Q c (y)dx 



(2) Almost solution. The error associated to the new function u(t) satisfies 
S[u] = (c'(t) - Eh(t))d c u - (p'(t) - c(t) - ef 2 {t))u y + S[u](t), 

with 

\\S[u](t)\\m m <e 3/2 e-^^. 

Finally, one has \3.11}) . 



(3.42) 
(3.43) 



(3.44) 



Proof. The proof of (|3.42[) follows from a direct computation. Indeed, 

\\w{t)n £ \\ H i {VL) < \\w(t)\\ H i(y>-Z), 

but from ([3~4"|) . 

\\ed{t)A c {y)\\m (y >^)<^e-^\. 
Let us now consider p.43|) . Here we have, using 



w(t,x)r) s (y)Q c (y)= / w(t,x)( V (ey + 2) - l)Q c {y). 
Js. 

Note that rj(ey + 2) — 1 = for y > — -. Using the exponential decay of Q c (y), we have 



w{t,x)r) £ {y)Q c {y) 



< 



e\y\e 



Vcy 



y<-i 

< e -7/e < 10_ 



| y | e -i(=y+2) e v^y 



The proof for yA c is very similar. We skip the details. For the proof of (|3.44p . we proceed as 
follows. First of all, a simple computation shows that 

S[r) s u] = i] £ S[u} + {n £ ) t u + 3sri' s u xx + "ie 2 r] [2) u x + e 3 ^ 3) u + e{rg)'vP. 

Since supp^^ C [— |, — i] for k = 1, 2 and 3, we have 

3e V ' e u xx + 3e 2 v Pu x + e 3 V ^u + e(rf E )'v? = O m{WL) {e 3 ' 2 e-^^) + H1{s) (e 10 ). 

Similarly, from the definition of p'(t) and (|2.5j) 

( Ve ) t u = - P '(t)e V ' £ u = Hm) {e*' 2 e-^rW) + O m{R) (e w ). 

Collecting the above terms, we have 

^[^ul^^^^+OHW^^e-^^^+O^^^ 10 ). 



Finally, from the decomposition (|3.8[) , one has S[u] = dynamical system + S[u], and from (|3.44j) , 
(pOl) and (f3~4l) . 

heS{u]\\H Hm <s 3/2 e-^+e 3 . 
Note that the proof of (|3 . 11 1) does not vary at all. Finally, one has 

r\ £ (c - efi)d c u + (p - c - ef 2 )d p u = (c - efi)d c (n £ u) - (p - c - ef2)d y (n £ u) 

+ e(p' - c- ef 2 )r\ £ u. 

Since en' £ u = O ff i( R )(£ 3 / 2 e _E7 ' p '*"), from this last estimate, we get the final conclusion. □ 
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Coming back to the original equation (|1.7[) - (|2.6[) . we have 

u t + {u xx + u v ) x + ea' (ex)Q Co (x - p )u = S[u] + ea' Q (ex)[Q Co (x - p ) - Q c {y)}u. 
With a little abuse of notation, we redefine 

S[u] := S[u] + ea' Q (ex)[Q Co (x - p ) - Q c (y)}u. 
Finally, using (|3.7|) . we get 

\\S[u](t)\\Hi W Ze 3/2 e- e7im . (3.45) 

4. Lyapunov stability 

So far, thanks to Propositions 12.11 and 13.71 we have well-defined the maps 
(co(i), Po(t)) satisfying ()2.5[) i — > u(t) mild solution of fj 1 . 1[) . 

and 

(c(t),p(t)) satisfying (|2.5p - (|3.2l) i — > u(t , c(t) , p(t)) approximate solution of ([l.ip , 

in the interval [0, T]. Now we will find a unique set of parameters (c, p), defined in the same 
interval, such that the corresponding functions u(t) and u(t) are close enough. In this section we 
prove the following 

Proposition 4.1. The following holds for any < e < Eq. There exist Kq > independent of e 
and unique C 1 functions c, p : [0, T] — > R such that, for all t <G [0, T], 

\\u(t)-u(t,c(t),p(t))\\ HHR) <K ^, (4.1) 

and 

\p'(t) - c(t) - ef 2 (t)\ + ±=\d{t) - eh{t)\ < KoV^- (4.2) 

Finally, one has 

| c (0)-l| + |p(0)+T|< Vi, (4.3) 

with constants independent of Kq . 

A direct conclusion of the previous result is the following 
Corollary 4.2. Let Tq > be the maximal time of existence of u(t). Then 

To > s' 1 ' 50 . 

Remark 4.1 (Notation). For the sake of brevity, in the forthcoming computations, we will denote 

c[ := c - sfi, and p[ := p' - c - e/ 2 . 

Proof of Proposition [4TT1 Let K* > 1 be a constant to be fixed later. Since u(0) — Q(x), by 
the local continuity in iJ 1 (M) of the gKdV flow, there exists a time T* > such that, for all 
t £ [0, T*], we can find continuous functions \(t),r(t) £ R, such that 

|K*)-«(- ,X(t),r(t))\\ Hlm <K*^, (4.4) 
\X(t) - co W| + \r(t) - P0 (t)\ < \\u(t) - u(. , A(t), r(t))\\^ m . (4.5) 

and 

|A(0)-l| + |r(0)| < Vi, 

with a constant independent of K* large. Without loss of generality, we can assume T* as the 
least upper bound of times such that the above properties are satisfied. The objective is to prove 
that we can take T* > T for K* large enough, by proving a bootstrap estimate for suitable well 
chosen parameters X(t), r(t). Our first step is to choose such parameters. 
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Lemma 4.3. Assume < e < Eq(K*) small enough. There exist unique C 1 functions c(t),p(t) 
such that, for all t G [0, T*], 

(4.6) 



(4.7) 
(4.8) 

(4.9) 
(4.10) 
(4.11) 



z(t) = u(t) — u(t, e(t), p(t)) satisfies / z(t, x)yQ c {y)dx — / z{t, x)Q c (y)dx = 0. 

Jn Jr 

Moreover, we have, for all t G [0, T*], 

||z(0)|| ff i (R) + |c(0) - 1| + |/3(0) + T| < Vi, 

ll*(*)l|in(R) + |c(*) - 03 (t) | + \p{t) - Po (t)\ < K*^fe. 
In addition, z(t) satisfies the following equation 

zt + {zxx + {u + z) m - u m } x + ea' (ex)Q ca z + S[u) + c[(t)d c u - p\{t)d y u = 0. 

Finally, there exists 7 > independent of K* such that for every t S [0, T*\, 

1/2 



\p[(t)\ 



< 



z 2 e - 7 ^|y| 



+ 



z\t) 



yQcS[u] 



and 



3 -7V5|»U2 



(n 



1/2 



Q c 5[u] 



Proof. The proof of (|4.6[) - (|4.7[) is a standard consequence of the Implicit Function Theorem, applied 
for each time t £ [0, T*]. Indeed, for fixed t G [0,T], let us define the map 

(v,c,p) G ff^M) x K + x R — >• J{v,c,p) G K 2 , 

where J = (Ji, J2) and 

(v-u(t,c,p))(x- p)Q c (x- p)dx, J 2 (v,c,p):= (v - u(t,c, p))Q c (x - p)dx, 



Jx(v,c,p) :-- 
It is clear that 

Ji(u(t, c (£),/9o(i))) c o(0,Po(i)) = J 2 {u(t 1 co(t),po(t)) 1 c (t),po(t)) = 0. 

Moreover, the respective Jacobian determinant of J with respect to the variables (c, p) is nonzero 
everywhere. Therefore, from the Implicit Function Theorem, there exists a small 770 (t) > (which 
can be chosen continuous on t) such that, for all v G H 1 (M.) satisfying \\v — u(t, Co(t), po (*)) II J? 1 (K) < 
770, there is a smooth pair of parameters (c(v), p(v)) G M 2 , satisfying J(v, c{v), p{v)) = 0. Since the 
interval [0,T*] is compact, we can ensure rjo > independent of t. 



Note that from (|4.4[) - (|4.5|) . the function u(t) satisfies, for each t G [0,T], 
\\u(t) - u(t,c (t),p (t))\\ H i {R) < K Q ^/e < 770, 



(4.12) 



provided e = £0(^0) is chosen even smaller. Therefore there exists a smooth pair of parameters 
(c(u(t)),p(u(t))) =: (c(t),p(t)) G R 2 , such that J(u(t), c(t), p{t)) = 0. This proves 



The proof of (|4.T[) is direct from the initial condition uo(x) = Q(x), and (|2.5p . Finally, 
follows from (|4. 12[) and the fact that J (u(t) , c(t) , p(t)) = 0. 



On the other hand, 



is a direct computation. For the proof of (|4. 1 0[) and (|4.1ip . see e.g. 

□ 



From (|3.1ip and (|4.7p . a crude estimate of the parameter |c^(t)| gives |c^_ (t) | < e, so that 



which is not a good estimate. In the following lines, we prove a sharp Virial estimate HZ] 
which allows to improve the preceding bound. 
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First of all, wc define some auxiliary functions. Let <fi € C°°(R) be an even function satisfying 
the following properties 



J (/)' < on [0, +oo); <j){x) = 1 on [0, 1], 

[ <f){x) = e~ x on [2, +oo) and e~ x < 4>{x) < 3e~ x on [0,+oo). 
Now, set ip(x) := Jo 4>- It is clear that ip an odd function. Moreover, for \x\ > 2, 

ip(+oo) - ^(\x\) = e"K 

Finally, for A > 0, denote 

^ A {x) := V(^); e-l^l/^ < ^ A (a) < 3e'^/ A . 

We claim the following 

Lemma 4.4. TTiere eirasf K, A 0l S > smc/i i/iai /or a/Z i € [0, T*], 

S t / z 2 {t lX )i>A {y) < S a I (z 2 x + z 2 )(t,x)er^ lyl +KA (K*r +1 e 5 / 2 . 



Proof. Let t G [0,T*]. Replacing the value of z t given by ()4.9|) . we have 
d t f z 2 ^ Ao (y)=2 [ zzti> A M-p'(t) I WaoM 



-(c + e/ 2 )(t) / z 2 iP' Ao -2p[(t) / z5 p «Va 

JR JR 

+2 / (z^o(y)),P + ^) p -« p -^ p "^] 

JR 

-2c£(t) f zd c ui/j Ao - p[(t) [ z 2 i/j' Ao 

JR JR 

+2p [ z^A^y))^- 1 - EP~ X ) - 2 f z^ Ao S[u] 

JR JR 

The first term in (|4.17j) can be estimated as follows: 

e / a' (ex)Q Co z 2 ip Ao (y) 



<e I z 2 (t)e~^ lvl . 



On the other hand, note that 



| (14~19D | < 



z x ipA (y) [(u + z) p -u p - pu p z 



(4.13) 

(4.14) 
(4.15) 



(4.16) 



2e / a' (ex)Q co z 2 ^ Ao (y)+2 \ (z^ Ao (y)) x (z xx + pBP~ 1 z) (4.17) 

(4.18) 
(4.19) 
(4.20) 
(4.21) 



^ Ao (y)z[(u + z)P~uP~puP- 1 z] 



SAo^e 1 / 2 I z 2 (t)e-~<^ +K*s 1 ' 2 I z 2 (t) e -^ lvl + 



<K*A e 1/2 I z 2 (t) e -^ M +A e\\z'^ p+1 



ff!(R) 



< K*A e^ 2 J z 2 (t) e -^ lvl + {K*)P+ 1 A e ( -P +3 '>/ 2 . 



for Aq large, but independent of e. Now, by using (|4.10|) and f)4. 1 1[) it is easy to check that for Aq 
large enough, and some constants Sq,£q small, one has 



f zd c iiip Ao +K*s 1 ' 2 f z 2 (t) 

JR JR 



< A. f z 2 (t)e-^ lyl + KK*A Q e 5 / 2 e-^rW. 
100 J K 
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On the other hand, the terms (|4.17|) and (|4.18l) goes similarly to the terms Bi and B 2 in Appendix 
B of [50]. Indeed, we have 



em + (En 



^ Aa (?zi+cz 2 -pQr i z i )-p / (or 1 )'* 3 ^ 



+ / z 2 ^] - 2p[(t) I zd p u^ Ao +2p I (^Mif-'-QT 1 )-^ / * 3 ^- 



We finally get, taking e small, depending on Aq, 



([4~T7]) + (j4~T8]) < 



So 
10 



(^ + « 2 )(t)e-^l*l. 



Finally, the term (|4.21[) can be estimated as follows 



< 



z(zipA (y))x(u P ~ 



R p - 



4>A S[U 

zz x tjjA {y){u p ~ 1 - i? p_1 ) 
< A e / (z 2 (t) + z 2 x (t))e'^ lyl + A K*s 5/2 e~ £ ~<^. 



< 



z^'aAvW- 1 - BP~ X ] 



+ A Q K*e^ 2 e-^M^ 



We have used that ipA decreases exponentially as y — >• — oo, and (|3.45[) . Collecting these estimates, 
we finally get (|3~Tg]) . □ 



Corollary 4.5. One /ias, /rom 1 1\ ) and 

I \c[{s)\ds<K*e, 
Jo 

for all t £ [0,T*], by taking Aq large enough, independent of e and K* . 



(4.22) 



The main part of the proof is the introduction of the following Lyapunov functional f |21[ 126] ): 



Let 



(4 + c(t)z 2 ) 



p+1 



\{u + z) p+1 - u p+1 - (p + l)u p z] 



(4.23) 



From [3T] and the fact that p < 5, there exists a constant, independent of K* and e such that for 
every t G [0,T*] 

Ht) > \\m\ 2 H H my (4-24) 
The next step is to obtain independent estimates on F(t). We follow [27]. It is not difficult to 
check that 



z t {z xx -cz + (u + z) p -u p ] + -c'{t) z 2 - u t [(u + z) p -u p -pu p - l z\. (4.25) 



Lemma 4.6. There exists a constant 7 > such that, for any t G [0, T*], 

ft 



T{t) - J"(0) < (iT)V~™5 + (x*) j £ 2-ttk7 + K*e + / ee 



-£7|p(«)l I 



z ( s )\\m(R)ds- 



Proof. Replacing ([43]) in (|4~2"5]) we get 

F\t) = e f a' (ex)Q Ca z[z xx - cz + (u + z) p - u p ] + c(t) I [(u + z) p - u p ]z x 



+c[(t) / d c u[z xx - cz + (u + z) p - u p ] 



+ / S[u][z xx -cz + (u + z) p -u p ] + -c' 1 (t) / z 2 + -ef 1 (t) / z 



iit[(u + z) p -u p -pu p - Y z\ 



(4.26) 

(4.27) 
(4.28) 
(4.29) 
(4.30) 
(4.31) 



IS 
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Wc consider the case p = 2, the other cases being similar (see [55] for more details). First of all, 
note that 



e / a Q (ex)Qc z[z xx - cz + (u + z) p - u p ] 
Jr 

Next, after some simplifications, we get 



<£e- 7£|p(t)l ||z(t)|' 2 



H 1 (R) ' 



c(t) I [(u + zf - u p ]z x = c(t) / [2uz + z 2 ]z x = -c(t) / u x z 2 . 
Jr Jr Jr 

Now we estimate (|4~28j) . Since d y u = Q' C + 0{r] £ w y ) + H i {m {e 3 / 2 e- £ ^P ( - t ^) (cf. Proposition EH]), 
one has 



Jr 



z xx — cz + 2uz + z 



-p'S) I u x z 2 + 0{ser^\\\z{t)\\l H 



(4.32) 

Jr 

Similarly, we have from (|4.6[) 

63S) = ci(t) ( d c u[z xx -cz + 2uz + z 2 ]=c\{t) [ d c uz 2 + 0(ee- £ ^\\z(t)\\ 2 L2m ). (4.33) 
Jr Jr 

On the one hand, we have from (|3.45p . 



S , [u]{z ;r:I , — cz + 2uz + z 2 } 



< 



d x S[u)z x 



+ (l + K*e 1 / 2 ) 









[ S[u]z 
Jr 


+ 


J S[u]uz 



< e 2 (e - £ 7|P(*)l + e )\\z(t)\\ m (M) + (1 + K*e^ 2 ) 
Concerning the second and third terms in (|4.30|) . 



S[u]z 



(4.34) 



c'i(t) / z 2 + -eh{t) / z 2 <(|c' 1 (i)l+« 



-£7|P(*)I 



)ll*(*)ll 



Finally, 



(fit + p'u x - c'd c u)z 2 + p' f u x z 2 -c' f d c uz 2 + 0{se-^Hz{t)\\ 2 L2m ), (4.35) 

Jr Jr 



$M = - 

We get then, 

^(*) < l^(*)llk(t)lli aW +^ 8 l"WI||z(t)||i a(R) + £ ||z(t)||^ (K) +^ £ 2 e-^l"«l. 
Using (|4.22p . we finally get after integration in time (here we use the condition t < T* < e^ 1 ^ 26 ") 

Hi) HO) < (in 3 e%- 5 ° + (K*s) 2 + [\e-^W\\z(s)\\% 1(m ds, 

Jo 

as desired. The cases p = 3 and 4 are similar. □ 



We are finally in position to improve (|4.4[) . Indeed, since from Lemma [4.31 J-(Q) < e, using 
(|4.24|) and Lemma (|4.26l) we get 

Mt)\\h m <e+(K*re 2 - So + (if*) 3 e^ 5 ° + {ICfe 2 + [\e-^^z( s )\\ 2 H1{R) d S . 

Jo 

Now, by Gronwall's inequality (see e.g. [25J for a detailed proof), 

\\m\\W)^ e +( K *) 3£i ~ So - ( 4 - 36 ) 

with constant independent of K* and e. Indeed, we just need to justify that J* ee~^ p ^ds < 1, 

independent of e and K* . It is not difficult to see that the above estimate holds since p'{s) > 
^c(s) > ^min{ C/ ,l} > 0. 
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Let us come back to the main proof. From estimate (|4.36[) . and taking e small, and K* large 
enough, we obtain that for all t G [0,T*], 

^ \( K *) 2 e- ( 4 - 37 ) 

Therefore, we improve the estimate on z(t) (|4.7[) , and therefore (|4.4p is also improved. The proof 
of Proposition |4~T1 is complete. □ 



5. The fixed point argument 

Let W be the nonlinear operator that associates the unique parameters (c, p) given in Lemma 
14.31 to the initial data parameters (c ,p ). We will prove that W has fixed point in a suitable 
small set of C 1 ([0,T],R + x R). Indeed, we consider the (convex) closed balls (recall that So > 
and < e < £o are small parameters) 

B := Uc°,p°) G C 1 ([0,r],R) : (c , p ) satisfy and 

sup |c°(t) - c (i)| + |p°(t) - PoWl < £ 1/2 ^ }, 
te[o,T] J 

and 

Bi == \(c,p) G CH^TJ.R) : sup |c(t) - co(*)| + W) - Po(t)\ < 

L te[o,r] J 

From Lemma T4.31 it is clear that W maps Bo into W^(i?o) C B± C _B . 

Proposition 5.1. The map W has a fixed point in Bq, provided Eo is small enough. 

Proof. It is not difficult to check that, thanks to Proposition 12.11 and Lemma 14.31 W defines 
a continuous map from Bo into itself. Let us prove that W(Bo) is a pre compact subset of 
Bo- Indeed, since (c, p) are diffcrentiable and from (|4.10j) - (|4.11j) one has W(Bo) is a uniformly 
bounded by 0(y/e). Moreover, W(B ) is clearly an cquicontinuous subset of B , endowed with 
the C 1 uniform norm. Thanks to the Schauder's Fixed Point Theorem, we conclude. □ 

So far, we have proved the existence of a couple of parameters (c,p), defined in [0, T], of class 
C 1 and such that Proposition 14. II is satisfied for the Cauchy problem (jTTTJ) with parameters (c, p). 
In the last section, we prove that such a solution satisfies the conclusions of Theorem 1 1.2 1 



6. Proof of the Main Theorems 



We are now in position to give a direct proof of Theorem 1 1.2 1 Indeed, we have (|4.ip for all time 
t G [0, T\. in particular, at t = T one has 

\\u{T)-u{T,c{T),p{T))\\ HKm < yfe. (6.1) 

Note that, from (|4~TT|) and (|4~T0|) . since c(t) > c by (|2~5|) . 

C -^(t)c'(t) = -e\ p a (ep{t))c(t) + {t)c[(t) 

= -e\ p a (ep{t))p'(t) + e 2 X p a' (ep{t))f 2 (t) + eA^M*))^*) +c~^{t)d l {t). 

Hence, if p = 2, and using f|4.22[) . 

log c(t) - logc(O) = -A 2 [ao(ep(t)) - ao(ep(0))] + 0(Vi), 

from which we obtain 

c(t) = e- A2a0 ^«)(l + O(V^)), p = 2. (6.2) 

Similarly, if p = 3 or 4, 

c(i) = [l - A p fc^a (ep(t))l ^ (1 + O(Vi)). (6.3) 
L p — 1 J 
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Now, we perform a detailed asymptotic analysis of (c(t),p(t)). First of all, note that from (|6.2[) - 
(16. 3[) . c(t) is strictly positive for all time. Indeed, the case p = 2 is direct, and for p = 3 or p = 4, 
we have c(t) > 1 — Ky/e in the case where doo is negative (or Cf > 1, sec (12.21) and (|2.3I) ). and 

JL 

for the case > (c/ < 1). From (|2.4[) we have ||ao||oo = floo = (1 — c| _1 ). Therefore, 

c(t) > Cf-Kyfe > 0. 

We conclude that, for e > small, c(i) > ^min{l,c/} > 0. Hence p(t) is increasing and 
p{t) - p(0) > ^ min{l,c which implies that p(T) > T. Moreover, from (E3J), a(ep(T)) = 
aoo + 0(e 10 ). Taking i = T in ([O]) and (JO|), we have 

c(T) = e" A2Q - (1 + O(Vi)) = c/(l + O(Vi)), p = 2, 

and 

■l 

3 (l + 0(Vi))=c/(l + 0(VS)), P = 3,4, 
as desired. Finally, from (|3.3[) , one has 

||«(T, c(T), p(T)) - Q c/ (• - p(T))\\ mm < y/i. 

Using (|6.ip and the triangle inequality, the first estimate in (| 1 . 8[) follows. Concerning the second 
one, it is a consequence of (|4.10p . Corollary 11.31 is just a consequence of the behavior of c(t) in 
([Pjl and 



=(*)> 



1 \ IP-^ n i 
1 - A p — — 1| cio I 



c(T) 



i — A„ — a a 

r> — \ 
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